Kinetic coarsening in the random field Ising model : a different approach 
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We report a study of nonequilibrium relaxation in a two-dimensional random field Ising model at 
a nonzero temperature. We attempt to observe the coarsening process from a different perspective 
with a particular focus on three dynamical quantities that characterize the kinetic coarsening. We 
give a simple generalised scaling theory of coarsening supported by numerical results. The excellent 
data collapse of the dynamical quantities justifies our proposition. The scaling relation turns out to 
be universal in characterizing the kinetic coarsening. 
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Study of the effect of disorder and randomness on 
nondisordered magnetic systems has been a subject of 
intense interest for the last several years [IHl]- When a 
system is quenched from a high temperature to a low tem- 
perature, it locally orders with the formation of domains 
separated by domain walls. The average linear size of 
the domains R{t) grows with time. This linear size can 
also be understood as the non equilibrium correlation 
length of the system. The growth of the characteristic 
length scale R{t) with time is known as the coarsening 
of the system. Although coarsening process in nondisor- 
dered systems is very well understood [loj . progress in 
understanding the same in disordered systems has been 
rather slow. In disordered systems, most aspect of non 
equilibrium relaxation remains under debate Un- 
availability of reliable theoretical tools makes it difficult 
to study the dynamics of disordered systems out of equi- 
librium. Moreover, the dynamics of disordered systems is 
typically so slow that we cannot access the truly asymp- 
totic time regime in numerical simulations. Despite all 
these, last several years have witnessed appreciable ef- 
fort in the study of disordered but unfrustrated systems. 
These include coarsening of disordered magnets jl3 - [l4l | , 
polymers in random media [T5l - [l7| or vortex lines in dis- 
ordered type - II superconductors [T8l-l20| . 

The ideal model for studying coarsening in disordered 
systems is the random field Ising model (RFIM) in which 
the random fields or impurities couple directly to the lo- 
cal order parameter. The fundamental quantity of in- 
terest in the coarsening process is the growing length 
scale R{t) and almost all studies of coarsening is pri- 
marily concerned with the determination of this charac- 
teristic length scale R{t). But the growth law govern- 
ing the coarsening process of disordered systems is at 
the center of some controversies. Some numerical sim- 
ulations on disordered ferromagnets (2ll - [23| yielded an 



algebraic growth R{t) ~ t^^^, with a nonuniversal dy- 
namical exponent z that depends on the tem per ature 
and on the nature of disorder. Huse and Henley [23| sug- 
gested a logarithmic increase of this characteristic length, 
R{t) ^ (Int)^/"^, with the barrier exponent tp > 0. Later 
a series of papers on the dynamics of elastic lines in a 
random potential [H, [13, HI] claims a dynamic crossover 
from a pre-asymptotic algebraic regime to a asymptotic 
slow logarithmic regime. 

In the present work, we attempt to study coarsening 
in a disordered spin system, namely the two-dimensional 
RFIM, from a different angle. To characterize the coars- 
ening, we focus our attention on three dynamical quan- 
tities which are, in general, functions of the strength of 
the random fields (770) and the temperature (T). These 
arc the total length of the interfaces {Il{ijo,t)), i.e., the 
total number of boundary spins of all the domains, the 
total number of domains (A(77o, t)) in the system and the 
length of the interface of the domain with largest mass 
(ri(77o,t)), i.e., the number of boundary spins of the do- 
main containing maximum number of spins. In general, 
we denote by "^{rjo,!) the three quantities which char- 
acterize the kinetic coarsening in the two-dimensional 
RFIM. 

The Hamiltonian of the RFIM is given by 
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where Si = ±1 is the spin variable at site z, J is the 
strength of the exchange interaction (conventionally set 
to unity) and r]i is the quenched random fields taken from 
an uniform distribution with varying strength 770. 
indicates that the sum is over nearest neighbours. The 
external field i?ext has been set to zero to observe the 
unbiased dynamics of the system. We consider a L x L 
square lattice (here L = 256) with periodic boundary 
conditions along both directions. We start our simu- 
lations from a completely random spin configurations, 
characteristic of a high temperature (T cx^) phase 
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and then suddenly queneh the system to a temperature 
T = 0.50, which is weh below the critical temperature 
of a nondisordered system (Ising model) and then ob- 
serve the time evolution of the system. The single spin 
flip Metropolis algorithm [2^ is used to simulate the sys- 
tem. Being a single spin flip dynamics, the Metropolis 
algorithm is believed to represent the natural way of the 
evolution of a system, since the acceptance ratio is given 
by the Boltzmann probability and the dynamics is lo- 
cal. The temperature is taken sufficiently low to reduce 
the thermal fluctuations. Coarsening proceeds through 
a compromise between the strength of the exchange in- 
teraction and the random fields with the thermal fluctu- 
ations ser ving only to renormalize the strengths of these 
couplings [27[. The number of domains with their sizes 
are determined by the Hoshen-Kopclman algorithm [28j . 
All the quantities arc averaged over 50 independent sim- 
ulations to get a precise estimate. The quantities are nor- 
malized with respect to the total number of spins (i^) of 
the system. 

We begin our analysis with the idea put forward by 
Imry and Ma In their introductory article, they ar- 
gued that if one reverses the spins within a domain of 
linear size i?, the energy cost E^-^ is proportional to the 
domain wall area, i.e., i?ex o^ JR'^^^ where d is the spa- 
tial dimension. This energy increase has to be compared 
with the energy gain from the interaction with the ran- 
dom field. The central limit theorem tells that the mean 
square random field energy i?^p inside a region of volume 
is ^-qlR^. The total energy involved in the creation 
of a domain of linear size R is therefore 
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The first term of Eq. ^ represents the contribution due 
to the boundary of the domain with linear size R(t). The 
second term represents the contribution due to the fluc- 
tuations of the random fields in the bulk of the domain 
of linear size R{t). We have extended this idea for deter- 
mining the energy corresponding to the creation of all the 
domains at any instant of time and for a particular disor- 
der strength (770). Taking a cue from the above relation, 
the surface energy of all the domains E^^ ^ jn(?7o, t) and 
the mean square bulk energy contained in all the domains 
of the system due to random fields E^^'^ ^ A{riot) ' 
the density of domains is inversely related to their char- 
acteristic volume. Thus the energy density (e = E^ /U^) 
involved in the creation of all the domains in the system 
is given by 



jn(77o,^) 
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As L ^ 00, the surface energy contribution vanishes. 
This is true for any growing volume. For finite system 
size, the contribution from the surface energy term can- 
not be neglected. So the variation of both n(?/o,0 and 
A(?7o,i)^^^ with time will govern the coarsening of the 



overall system. The log-log plot of Il{T]o,t) and A{r]o,t) 
against time are shown in Fig. [l]and Fig. [5] respectively. 
It is evident from Eq. ([3]) that the coarsening of the 




FIG. 1. The plot of the function 11(770, against time. The 
dotted lines are the best fits according to the scaling relation 
given by Eq. ((?]). 




FIG. 2. The plot of the function A{rio,ty^' against time. The 
dotted lines are the best fits according to the scaling relation 
0. 

system energetically favours the minimization of the to- 
tal length of the interfaces and also the decrease of the 
number of the domains. This is observed in FiglT] and 
Figl2] respectively. In view of the above discussions, we 
can redefine the coarsening as simply the minimization of 
e(riQ,t) and the kinetic coarsening will be characterized 
by the scaling behavior of the two quantities 11(770, t) and 
A(77o, ty^^. With the fiow of time, small domains coalesce 
to form relatively larger domains and from Eq. it is 
clear that a domain with a typical size R should grow in 
such a way that the length of the interface of that domain 
shrinks while the volume of the the domain expands in 
order to minimize the energy of the domain. The domain 



with largest mass should grow in the same fashion during 
its dynamical evolution and plays the most dominating 
role in governing the dynamic behavior of the system. 
Recently, we report in [29| that the size of a domain in 
two-dimension increases as a power law with time. Since 
the domain with largest mass predominates over all the 
other domains, the length of the interface of this domain 
{Cl{rio, t)) is expected to exhibit similar behaviour as that 
of n(?7o, t). Log-log plot of il{T]Q,t) against time is shown 
in Fig. H 
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FIG. 3. Plot of Q,{'qo,t) against time along with the best fits 
according to ((T)). 

We now give a simple unified coarsening theory. A 
careful observation of the graphs suggests that the ini- 
tial and asymptotic behaviour of the generalized function 
^{rjQ,t) is given by 

^(770, — > ^"0 (constant), as t — ^ 1 

and ^'(?7o,i) -J- ^'oe"^/''(''°\ i'(7?o) > 0, as t ^ o^) 

where i'(f?o) is a disorder dependent scaling exponent. 
The decay rate of the function »I'(7yo,i) at any time step 
for a particular value of rjQ should depend on the following 
factors: first, on the value of the function itself at this 
time step. Secondly, from the nature of the variation 
of the functions Il{r]o,t), A{r]o,ty^'^ and 51(770, i), it is 
evident that the rate of decay of the function \E'(r/o,i) 
also depends on the particular time step. As the time 
flows, the rate of decay of 'if{r]Q,t) slows down and this 
dependence is taken as a power law decay. In addition 
to these factors, another ryo-dcpendcnt parameter should 
come into the picture for controlling the decay rate of 
\E'(77o,t). This parameter considers the wandering of the 
interfaces in presence of the random field. Thus the decay 
rate of 4'(r?o, is given by 

— ^ -a(77o)*t-''<*\ ^^M>0 (5) 
at 

0(770) is a disorder-dependent parameter. Integrating, 

"f-(M(»?o)-l) 



«'(7/o,t) = foexp 
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where A; (770) is a constant of integration and 7^(770) = 
tin^yiL Now from © as t 1, ^(770, *o which 
1/1^(7^0) and as i — > 00, 'i/{ijo,t) — > 



gives fc(77o) 
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which gives ^{rjo) > 1. Thus the functional 



form of 5'(77o,t) at any time step is given by 



^-(770,0 = *oexp 
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where p(7;o) = /z(7;o) — 1 > 0. The scaling behaviour (O of 
the functions characterizing the kinetic coarsening shows 
an universal nature with two disorder-dependent expo- 
nents p(77o) and v{r]o). The plots of the data collapse for 
the three functions 11(770, i), A(77o,t) and 11(770, i) charac- 
terizing the kinetic coarsening are shown in Fig. 21 Fig. 
[5] and Fig. [H] respectively. 
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FIG. 4. Plot of the data collapse of the function 11(770, t). The 
inset shows the variation of p(77() and !^(»7o) against 770. 
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FIG. 5. Plot of the data collapse of the function A{rio,ty^^ . 
The inset shows the variation of p(77o and 1^(770) against 770. 

It is evident from these figures that the scaling rela- 
tion (O yields an excellent data collapse, which confirms 
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FIG. 6. Plot of the data collapse of the function Q,{r]o,t). The 
inset shows the variation of p(r;o and !^(r;o) against rjo- 

the validity of our proposition. We interpret 1^(770) to 
be comprised of two parts. One part (say i^o) is fixed 
and it is a measure of tlie surface tension wlien the do- 
mains grow due to exchange interaction. The second 
part (say z^i(?7o)), that depends on i]q prevents the do- 
mains to grow with time. This part may be thought to 
be a measure of the stiffness of the interfaces. Beyond a 
threshold value of 770, i^iirjo) predominates after the time 
scale tx ~ l/(z/(?7o))^/''^'"'-' and the domains are pinned, 
i.e., they are prevented from growing. For t » t^, the 
dynamic function ^'(770,^) becomes almost constant with 



time and the asymptotic value is ~ 'i'oe~^^'^^'^°\ With 
this interpretation, 7^(779) should be a monotonically in- 
creasing function of 770- It is shown in the inset of the 
figures. The variation of p{rio) against 770 is also shown 
in the inset of the figures. As t — > 1, 5'(77o,<) reaches a 
fixed value 4*0, independent of 770. 

In summary, we have studied the kinetic coarsening 
phenomena in a disordered but unfrustrated system, 
namely the RFIM in a different approach. The impor- 
tance of this study lies in the fact that disorders are inher- 
ent in many physical systems and therefore the behaviour 
of these systems far from equilibrium is an active area of 
basic research in its own right. Most of the conventional 
studies of coarsening focuses considerable attention on 
the growth of the average linear domain size following 
a quench from a high temperature to a very low tem- 
perature. In the present work, instead we have devoted 
most attention on the time evolution of three dynamical 
quantities, introduced earlier in this paper. These three 
quantities characterize the kinetic coarsening of the sys- 
tem. We propose a dynamical scaling relation for the 
three quantities and have shown that the scaling is uni- 
versal in characterizing the kinetic coarsening of the two- 
dimensional RFIM. We hope that the approach to study 
coarsening along the lines done here for the RFIM can 
also be followed in other disordered systems as well and 
the same dynamical scaling holds good. 
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